Isospin-tracing: A probe of non-equilibrium in central heavy-ion collisions 
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Four different combinations of ftRu and foZr nuclei, both 
as projectile and target, were investigated at the same bom- 
barding energy of 400A MeV using a 4tt detector. The degree 
of isospin mixing between projectile and target nucleons is 
mapped across a large portion of the phase space using two 
different isospin-tracer observables, the number of measured 
protons and the t/ 3 He yield ratio. The experimental results 
show that the global equilibrium is not reached even in the 
most central collisions. Quantitative measures of stopping 
and mixing are extracted from the data. They are found to 
exhibit a quite strong sensitivity to the in-medium (n,n) cross 
section used in microscopic transport calculations. 

PACS numbers: 25.75.-q; 25.75.Dw; 25.75.Ld 

Central heavy-ion collisions represent a valuable tool 
for studies of hot and dense nuclear matter where one 
hopes to infer valuable information on the nuclear equa- 
tion of state (EOS) and on modifications of hadrons in the 
nuclear medium. R is still an open question whether the 
widely applied, at least local if not global, equilibrium 
assumption is valid in such reactions [[[J@], or whether 
significant non-equilibrium effects rather require the ap- 
plication of more elaborated non-equilibrium dynamical 
models [§-j4|- The issue of equilibration is expected to 
be influenced by in-medium effects (such as Pauli block- 
ing, Fermi motion) on the 'hard' scattering processes, by 
early 'soft' deflections in the momentum-dependent mean 
fields, and by finite-size (corona) effects. An understand- 
ing of all these effects is a prerequisite for a quantitative 



extraction of the EOS from nucleus-nucleus collisions. 

Experimental observations of non-equilibrium in rela- 
tivistic heavy-ion collisions were concentrated, up to now, 
on the measurement of the momentum distribution of the 
products emerging from a mid-rapidity " source" of sym- 
metric colliding systems. Observables of interest were the 
width of rapidity distributions j^J^] or the overall shape 
of the source The sensitivity of such observables 

is however reduced by effects like rescattering during the 
late phase of expansion. 

In order to extract, in a model independent approach, 
direct experimental information on non-equilibrium we 
have designed a new type of high precision measurement 
which makes use of the isospin (N/Z) degree of freedom. 
The (N/Z) equilibration has been investigated before at 
low bombarding energies in fusion- like reactions [p[-jIo[. 
As it will be shown in some details later, (N/Z) can be 
used as a tracer in order to attribute the measured nu- 
cleons (on average) either to the target or to the projec- 
tile nucleons. It is therefore possible to extract rapidity- 
density distributions separately for projectile and target 
nucleons. This gives access to new more sensitive observ- 
ables, like stopping and mixing, of the early equilibration 
process. 

The experiment was carried-out using reactions be- 
tween equal mass nuclei A = 96, at an incident energy 
of 400A MeV. Isotopes of Ru and Zr were taken as pro- 
jectile and target making use of all the four combina- 
tions: Ru+Ru, Ru+Zr, Zr+Ru and Zr+Zr. This choice 
of isotopes takes advantage of an almost unique possibil- 
ity offered by the periodic table of stable isotopes, while 
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searching for two isobars of the largest possible N/Z dif- 
ference (the N/Z ratio is equal to 1.18 and 1.40 for 44RU 
and 40 Zr, respectively) which can be used both as pro- 
jectile and target. The bombarding energy was chosen at 
the minimum of the (n,n) cross section where the relative 
motion is significantly larger than the Fermi motion, but 
sufficiently low to avoid inelastic (n,n) channels, while 
the (n,n) angular distribution is almost isotropic. 

The experiment was performed at SIS/ESR-Darmstadt 
using the FOPI apparatus jll],[l2| ■ The reaction products 
were detected at forward angles by the TOF-wall and the 
HELITRON drift chamber and at backward angles by the 
Central Drift Chamber (CDC). Particles of interest here 
were: i) tritons and 3 He emitted at laboratory angles 
10° < 9 lab < 28°, identified by combining the TOF-wall 
and the HELITRON and ii) protons and deuterons emit- 
ted at backward angles (34° < #i a b < 145°), identified 
by the CDC. These particles constitute a large fraction 
of the emitted charge, about 40%, 30%, 10% and 10% for 
p, d, t and 3 He, respectively. 




FIG. 1. Invariant cross section in the (y , p\ ) plane 
of Z = 1 particles measured in the CDC and the TOF-wall, 
under the -Erat > 1.4 centrality cut. j/ - 1 is the normalized 
rapidity, i.e. the particle cm. rapidity divided by the projec- 
tile rapidity in the cm. system. p[ ' denotes the normalized 
transverse momentum, i.e. the particle transverse momen- 
tum per nucleon divided by the cm. projectile momentum 
per nucleon. The grey levels correspond to different loga- 
rithmic cuts in the invariant cross section. The dashed lines 
show the acceptance of the CDC (34° < (9i a b < 145°) and 
the HELITRON (10° < 6 uh < 26°) detector components. 
The circles correspond to constant cm. energies per nucleon: 
from the inner to the outer, 40, 80, 120, 160 and 200A MeV. 

Central collisions were selected requiring a large ra- 
tio of total transverse to total longitudinal energy, -Erat 
HQ. Erat is determined event wise by including all 
detected particles: Erat = Si=i ^t/Si=i wnere E\ 
and E\ are the transverse and longitudinal kinetic ener- 
gies of the particle i and M is the number of detected 
particles. The Erat > 1.4 cut used here corresponds to 
about 1.5% of the total reaction cross section and to a 
geometrical impact parameter 6 goom < 1.3 fm in a sharp- 



cut-off approximation. 

Fig. 1 shows that the events selected under the 
Erat > 1-4 condition exhibit, in a representation of ra- 
pidity versus transverse momentum, a nearly isotropic 
source for the observed hydrogen (Z = 1) products. This 
is in agreement with previous observations [p|-|7|]l3|] . 

The (iV/Z)-tracer method is based on the following 
idea: let us assume that we are observing the final num- 
ber of protons, Z in a given cell of the momentum space. 
The expected yield Z Ru measured for the Ru+Ru reac- 
tion is higher than Z Zr of the Zr+Zr reaction since Ru has 
44 protons as opposed to 40 for Zr. Such measurements 
using identical projectile and target deliver calibration 
values Z Ru and Z Zr for each observed cell. In the case of 
a mixed reaction, Ru+Zr or Zr+Ru, the measured proton 
yield Z takes values intermediate between the calibration 
values (Z Ru , Z Zr ). If e.g. Z is close to Z Ru in a Ru+Zr 
reaction, means that the cell is populated predominantly 
from nucleons of the Ru-projectile while if it is close to 
Z Zr it is mostly populated from nucleons of the Zr-target. 
In this way it is possible to trace back the relative abun- 
dance of target to projectile nucleons contributing to a 
given cell. If the mixed value Z varies linearly between 
the two extremes of the calibration values (Z Ru , Z Zl ) as 
a function of the relative abundance of target versus pro- 
jectile nucleons, such an abundance can be most easily 
derived. As it will be shown this is indeed the case for 
protons but it is not the case if we use as a (JV/Z)-tracer 
variable not the number of protons but e.g. the relative 
tritium to 3 He abundance in the cell. The t/ 3 He-ratio 
varies non-linearly between the two extreme of the cali- 
bration reactions, but this can be taken into account and, 
as it will be shown, consistent results are obtained with 
different (iV/Z)-tracer variables. 

The following definition for the relative abundancy of 
the projectile-target nucleons has been adopted: 

2 x Z - Z Zr - Z Ru 



where Rz takes +1 for Zr+Zr and —1 for Ru+Ru. In the 
case of full equilibrium in a mixed reaction, Rz would be 
equal to everywhere independently of the location of 
the cell (for the case of a linear dependence of Z) . 

Several important advantages of the method can be 
mentioned: (i) the four reaction combinations are inves- 
tigated, under identical experimental conditions so that 
the ratios are insensitive to systematic uncertainties due 
to the apparatus. The errors are essentially of statis- 
tical nature and profit from the high yield in the con- 
sidered cell, (ii) the mixed reaction Zr+Ru is the same 
as Ru+Zr except that target and projectile are inverted: 
this allows forward-backward cross-checks of the appara- 
tus which in addition can also be obtained from the sym- 
metric Ru+Ru and Zr+Zr reactions, (iii) using the four 
reactions the full information needed can be obtained by 
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measuring only within the center-of-mass (cm.) back- 
ward or only the c.m. forward hemisphere. 
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FIG. 2. Two isospin-tracer observables as a function of 
the geometrical impact parameter: Rz (panel a) and R t pn c 
(panel b). They are determined in the c.m. rapidity range 
-0.75 < y {0) < -0.25 and 0.25 < y (0) < 0.75, respec- 
tively. The horizontal dotted lines correspond to the values 
(Rz = 0.04 and -R t/ /3 Ho = —0.27) at halfway between the 
Ru+Zr and Zr+Ru experimental points in the 6 gC om < 1.5 fm 
region. The dotted-dashed lines are just to guide the eye. 

The results to be presented here concentrate on two 
tracer-observables: (A) the sum (Z) of the number of de- 
tected free protons plus the number of protons detected 
within the deuterons. (B) the relative abundance of tri- 
tium to 3 He (t/ 3 Hc). Observable (A) was measured in 
the backward c.m. hemisphere using the CDC: this has 
the advantage to minimize background originating from 
the oxygen content of the (Zr02) target. Observable (B) 
was measured in the forward c.m. hemisphere using HE- 
LITRON and TOF-wall profiting from a good time-of- 
flight particle identification. 

Results on the centrality dependence of the abundancy 
ratios are presented in Fig. 2. The selected momen- 
tum cell is rather wide and integrates over a large rapid- 
ity bin. Panels (a) and (b) show results obtained from, 
proton yields and t/ 3 He-ratios, respectively. The quan- 
tity i? t /3 Hc is defined in a similar way as for Rz, using 
equation (1) with Z being replaced by the t/ 3 He abun- 
dancy ratio. The impact parameter 6 gcom is derived by 
integrating over the measured cross section as a function 
of .Erat or of the charged particle multiplicity. Except 
for an off-set of R t/ 3 Bc = -0.27 ± 0.07 for the t/ 3 Hc- 
ratio, both figures display the same trends. Such an off- 
set can be understood from a non-linear dependence of 
the t/ 3 He ratio as a function of isospin. We find that, 
for N > Z, an empirical dependence of the t/ 3 He ratio 
/ \ 2 ' 5 

of [1 + 100 x ( jjt^ ) ] describes the known system- 

atics (see Fig. 35 in ref. |14|) and satisfies the measured 
off-set of —0.27. The off-set on Fig. 2a is very small 
(0.04 ± 0.03) and demonstrates the linear dependence of 
variable (A), what is also confirmed in Fig. 3a: at mid- 
rapidity = 0) the ratio Rz is very close to 0, for 
both mixed reactions as expected from symmetry. 



The results of Fig. 2 are consistent with the picture 
of an incomplete global equilibrium which persists up to 
the end of the collision, an observation that holds even 



for the most central collisions, for b„ 
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FIG. 3. Panel (a) shows the tracer observable Rz as a 
function of the normalized c.m. rapidity for central collisions. 
The experimental results are shown for both isospin asymmet- 
ric reactions Ru+Zr (full circles) and Zr+Ru (open circles). 
The error bars correspond to statistical uncertainties. The 
horizontal dashed lines indicate the values of Rz = +1 and 
— 1 corresponding to the isospin symmetric reactions Zr+Zr 
and Ru+Ru, respectively. Panel (b) displays the experimental 
c.m. rapidity distribution (circles) of protons, free plus those 
detected within the deuterons, in central Ru+Ru collisions. 
Data points in the backward hemisphere (full circles) were ob- 
tained from the CDC detector. Those in the forward region 
(open circles) were obtained by assuming a backward/forward 
symmetry. The two other distributions have been obtained by 
unfolding the overall distribution into 'projectile' (squares) 
and 'target' (triangles) components. 

After the inspection of the centrality dependence of 
Fig. 2, we show in Fig. 3 the rapidity dependence of 
Rz for the most central collisions. For this purpose we 
use the proton yield (A) whose relatively high cross sec- 
tions permit a fine subdivision of the rapidity bins. In- 
verting projectile and target changes, as expected, the 
sign of the i?z-values; besides that, both results agree 
within errors. The dashed-dotted line (±0.393 ?/ 0) ) 
describes an average of both measurements (being the 
same except for the sign) and can be used to deconvo- 
lute the overall measured rapidity distribution diV/d?/ ) 
(Fig. 3b) for Ru+Ru into separated rapidity distribu- 
tions for the projectile- and for the target-nucleons. For 
each rapidity bin, the number of projectile (target) nu- 
cleons was obtained as iV pr = 0.5 (1 + 0.393 y (0) ) N 
(N tI = 0.5 (1 - 0.393 2/(°)) N). The overall dN/dy^ 
distribution was obtained by extrapolating the measured 
transverse momentum spectra, within the backward de- 
tector acceptance, according to the procedure described 
in ||. After deconvolution a shift between the two de- 
duced rapidity distributions emerges, demonstrating that 
a memory of the initial target /projectile translatory mo- 
tion survives throughout a central collision. 
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The separated projectile (or target) rapidity distribu- 
tion can be readily parametrized into: (i) a shift (yp? ) of 
the mean value with respect to mid-rapidity; (yp?) = 
corresponds to full stopping and full thermo/chemical 
equilibrium, while (ypr) = 1 corresponds to the ini- 
tial projectile rapidity without any stopping; positive 
values of (y$) are expected for transparency, nega- 
tive values for a backward rebound of the projectile 
nucleons from the target, (ii) a mixing value M pr = 
(Nj — Nb)/(Nf + Nb), where Nf is the number of pro- 
jectile nucleons emitted forward and Nb backwards; M pr 
is a measure of non-equilibrium effects, (iii) a width of 
the unfolded distribution, er pr . 

TABLE I. Mean value (ypr), width a pr of the projec- 
tile rapidity distribution (in ?/ ^ units) and mixing observ- 
able M pr (see text) obtained for Ru+Ru. Systematic errors 
on the values of (ypr) and M pr extracted from the data are 
about 10%. (1): IQMD predictions obtained by applying the 
(JV/Z)-tracer method. (2): same as (1) but using the de- 
tector filter. (3): same as (1) but here IQMD predictions 
are obtained by tagging projectile and target nucleons in the 
model. (4): same as (3) with a™ d = 0.5 crl^ c . (5): same as 
(3) with cC od = 1.2 <T^ e . (6): same as (3) but without MDI. 
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(n,n)-cross section, <r"J 1 cd , is shown in Table I. Both ob- 
servables (yp?) and M pr exhibit a quite strong sensitivity 
to <r™° d and depend also on MDI effects. On the other 
hand, the influence of the stiffness of the EOS was found 
to be very weak: calculations using a soft EOS (K = 200 
MeV) led, within statistical uncertainties, to the same 
results. A comparison to the experiment leads, within 
this model, to a value of crJJ I1 ed slightly higher (by about 
20%) than the free (n,n)-cross section, of'™: a significant 
reduction of cr I1 J 1 ed seems to be excluded. 

In conclusion, the proposed experiment and method 
demonstrate a high sensitivity to non-equilibrium effects, 
which are found to persist up to the end of a reaction 
even in the most central collisions. This confirms the 
necessity to use non-equilibrium dynamical calculations 
for studies of in-medium effects and of the EOS of high 
density nuclear matter. 
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The defined values are well suited in order to char- 
acterize with few numbers the separated distributions 
and the strength of the non-equilibrium effects. They 
are useful in order to verify with the help of theoretical 
models jis] the relevance of the (AT/.Z) -tracer method 
and effects due to the experimental filter. Such studies 
have been performed here in the context of the IQMD- 
model 15 1 . Calculations were done using a stiff (K = 380 
MeV) EOS and momentum dependent interactions (MDI) 
for impact parameter selections similar to those applied 
to the data. The results (Table I) show that the ex- 
tracted values are reliable. As shown from the IQMD 
calculations, the effect of the detector filter is negligible. 
This is due to the fact that the method relies on relative 
quantities. The systematic errors (~ 10%) on the experi- 
mental (j/p? ) and Mp r values take into account the influ- 
ence of other particle species (t,a and heavier fragments) 
not considered in the analysis. It is interesting to note 
also in Table I that the values obtained by applying the 
(iV/Z)-tracer method to the IQMD events are very close 
to those extracted by tagging the projectile and target 
nucleons in the model. This confirms the validity of the 
(iV/Z)-tracer method proposed in this Letter. The de- 
pendence of the extracted parameters on the in-medium 
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